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•  A  2-D  model  for  thermoelectric  cells  is  advanced  based  on  the  finite-volume  method. 

•  The  model  accounts  for  the  Fourier,  the  Thomson,  and  the  Joule  effects. 

•  The  model  predictions  were  validated  against  experimental  data. 

•  The  effects  of  the  thermoelectric  properties  and  the  cell  geometry  were  assessed. 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  12  March  2014 

Received  in  revised  form  20  May  2014 

Accepted  21  May  2014 

Available  online  13  June  2014 


Keywords: 

Thermoelectric  cooling 
Two-dimensional  model 
Sensitivity  analysis 
Finite-volume  method 


The  present  paper  is  aimed  at  putting  forward  a  two-dimensional  model  for  thermoelectric  cells.  The 
energy  conservation  equation  was  formulated  in  order  to  account  for  the  Fourier  (heat)  conduction, 
the  Thomson  (thermoelectric)  effect,  and  the  Joule  heating  on  the  temperature  distribution.  The  electric 
field  was  also  modelled  in  order  to  come  out  with  the  current  and  voltage  distributions.  The  governing 
equations  were  discretized  by  means  of  the  finite-volume  method,  whereas  the  TDMA  algorithm  was 
adopted  for  solving  the  sets  of  linear  equations.  An  explicit  solution  scheme  was  employed  to  address 
the  temperature  influence  on  the  thermoelectric  effect.  The  model  results  have  been  compared  with 
experimental  data,  when  a  satisfactory  agreement  was  achieved  for  both  cooling  capacity  and  COP,  with 
errors  within  a  10%  band.  In  addition,  the  model  was  employed  to  assess  the  effects  of  the  thermophysical 
properties  and  the  couple  geometry  on  the  thermodynamic  performance  of  the  thermoelectric  cell. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

In  the  past  decades,  solid-state  cooling  technologies  have  come 
onto  some  particular  market  niches,  especially  the  applications 
related  to  portable  cooling  [1].  The  most  significant  advances  have 
been  achieved  in  the  realm  of  thermoelectric  cooling  [2],  in  which 
an  electric  current  produces  a  temperature  difference  in  a  couple  of 
dissimilar  semiconductor  materials.  A  typical  thermoelectric 
module  is  manufactured  with  two  thin  ceramic  wafers  and  an 
array  of  p-  and  n-type  blocks  of  doped  semiconductor  material 
sandwiched  between  them.  A  pair  of  p-  and  n-type  blocks 
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connected  electrically  in  series  and  thermally  in  parallel  make  up 
a  thermoelectric  couple  [3]. 

Several  studies  have  been  conducted  both  theoretically  and 
numerically  to  assess  the  thermodynamic  performance  of 
thermoelectric  cells.  Samples  of  the  most  influencing  works  are 
summarised  in  Table  1.  The  literature  review  points  out  that  most 
models  are  one-dimensional,  being  not  able  to  evaluate  the 
influence  of  the  cell  geometry  on  its  performance.  In  addition, 
the  literature  analysis  also  reveals  that  the  few  available  multidi¬ 
mensional  approaches  are  often  developed  aided  by  commercial 
packages,  which  provide  restricted  access  to  the  mathematical 
formulation  and  the  numeric  algorithm.  At  last,  most  models  do 
not  account  for  the  heat  transfer  in  the  air  cavity,  which  may  also 
affect  the  system  performance.  The  present  paper  is  therefore 
aimed  at  advancing  a  tailor-made  two-dimensional  model,  in  the 
realm  of  the  irreversible  thermodynamics,  which  is  suitable  to 
evaluate  the  sensitivity  of  the  thermophysical  properties  and  the 
cell  geometry  on  its  thermodynamic  performance. 
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Nomenclature 

Z 

figure-of-merit  (K  *) 

Roman 

COP 

coefficient  of  performance  (WW-1) 

Greek 

j 

electric  current  density  (A  m-2) 

oc 

Seebeck  coefficient  (V  I<-1) 

k 

thermal  conductivity  (W  m-1  K-1) 

y 

electrical  conductivity  (A  V-1  m-1) 

Lx 

width  (m) 

(p 

generic  variable 

Ly 

height  (m) 

p 

electrical  resistivity  (Vm  A-1) 

Lz 

length  (m) 

T 

Thomson  coefficient  (V  I<-1) 

m 

number  of  control  volumes  (y-direction) 

w 

dimensionless  response  variable 

n 

number  of  control  volumes  (x-direction) 

2 

coefficients  of  Eq.  (27) 

N 

number  of  thermoelectric  couples  in  the  cell 

V 

dimensionless  values  of  (p 

Nu 

Nusselt  number  (dimensionless) 

q 

volumetric  heat  generation  rate  (W  m-3) 

Subscripts 

q 

heat  flux  (W  m-2) 

c 

cold  end 

Q 

heat  transfer  rate  (W) 

e,  w,  n,  s 

control  surfaces 

T 

temperature  (K) 

P,  E,  W,  N,  S  control  volumes 

V 

voltage  (V) 

h 

hot  end 

2.  Mathematical  model 

A  thermoelectric  cell  is  comprised  of  several  couples  of  p  and  n 
semiconductors  connected  electrically  in  series  and  thermally  in 
parallel,  and  separated  from  each  other  by  a  cavity  filled  with  air. 
In  the  present  work,  the  physical  model  is  restricted  to  a  thermo¬ 
electric  couple,  as  illustrated  in  Fig.  1,  which  in  turn  is  subdivided 
into  ten  domains,  as  summarised  in  Table  2.  The  dimensions  in 
Table  2  refer  to  the  thermoelectric  device  under  analysis  13], 
which  has  been  taken  as  reference  for  the  present  study  given  its 
application  to  portable  coolers  [14]. 

The  mathematical  model  is  based  on  the  following  key  assump¬ 
tions:  (i)  steady-state  two-dimensional  approach,  (ii)  thermoelec¬ 
tric  properties  regarded  as  functions  of  the  temperature  only,  (iii) 
negligible  internal  contact  resistances  (both  thermal  and  electric), 
(iv)  both  n  and  p  elements  have  the  same  Seebeck  coefficient,  but 
with  different  signs,  and  (v)  heat  transfer  by  both  advection  and 
radiation  are  disregarded,  so  that  Nu  =  1  in  the  cavity.  Thus,  a  local 
energy  balance  yields, 

Vq  =  q  (1) 

where  q  is  the  rate  of  heat  generation,  and  the  heat  flux,  q,  is  calcu¬ 
lated  from  the  following  fundamental  relation  obtained  from  the 
irreversible  thermodynamics  15]: 

q  =  -kVT  +  gcTj  (2) 

where  first  term  on  the  right-hand  side  stands  for  the  heat 
conduction  (referred  hereafter  as  Fourier  effect),  where  k  is  the 
thermal  conductivity,  and  the  second  term  is  associated  with  the 


+ 


Fig.  1.  Schematic  representation  of  the  physical  model. 

Seebeck  effect,  where  a  is  the  Seebeck  coefficient.  The  divergent 
of  Eq.  (2)  yields, 

V  •  q  =  -V  •  (fcVl)  +  Tj  ■  Va  +  a/-  V7  +  al(v  ■/)  (3) 

where  V  j  =  0  at  steady-state  conditions  to  ensure  the  continuity 
of  the  current  density.  In  addition,  the  definition  of  electric  field 
yields, 

-  W  =  pj+  ocVT  (4) 


Table  1 

Recent  literature  on  performance  assessment  of  thermoelectric  cells. 


Author 

Year 

Approach 

Thomson  effect 

Cavity  convection 

Cavity  radiation 

Physical  domain 

Properties  as  f=j[T) 

Huang  et  al.  [4] 

2005 

Analytical 

Yes 

Yes 

Yes 

ID 

No 

Pramanick  and  Das  [5] 

2006 

Analytical 

Yes 

No 

No 

ID 

No 

Lee  and  Kim  [6] 

2006 

Numerical 

No 

No 

No 

ID 

No 

Yamashita  [7] 

2009 

Analytical-experimental 

Yes 

No 

No 

ID 

Yes 

Chen  et  al.  [8] 

2011 

Numerical 

Yes 

Yes 

Yes 

3D 

No 

Meng  et  al.  [9] 

2011 

Numerical 

Yes 

Yes 

Yes 

ID 

Yes 

Du  and  Wen  [10] 

2011 

Numerical-experimental 

Yes 

No 

No 

ID 

Seebeck  only 

Chen  et  al.  [11] 

2012 

Numerical 

Yes 

No 

No 

3D 

Seebeck  only 

Perez-Aparicio  et  al.  [12] 

2012 

Numerical 

Yes 

Yes 

Yes 

3D 

Yes 
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Table  2 

Summary  of  the  physical  subdomains  and  thermophysical  properties. 


Subdomain 

Description 

Material 

Dimensions  (mm) 

k  (W/mK) 

y  (V/Am) 

1  and  6 

Bottom  (1)  and  top  (6)  electric  insulators 

A1203 

4.8  x  0.62 

30 

_ 

2  and  9 

Left  (2)  and  right  (9)  electric  conductors  at  the  bottom 

Cu 

1.9  x  0.41 

400 

5.96xl07 

3  and  7 

n-type  (3)  and  p-type  (7)  semiconductors 

Bi2Te3 

1.4  x  1.14 

Eq.  (12) 

Eq.  (11) 

4 

Electric  conductor  at  the  top 

Cu 

1.9  x  0.41 

400 

5.96  x  107 

5  and  10 

Left  (5)  and  right  (10)  side  air  cavities 

Air 

0.5  x  1.55 

0.026 

- 

8 

Central  air  cavity 

Air 

1.0  x  1.55 

0.026 

- 

Therefore,  the  rate  of  heat  generation  q  can  be  calculated  from: 

Q  =!■(-' Vv)=p(j-j)+ocjVT  (5) 

where  the  first  term  on  the  right-hand  side  stands  for  the  Joule 
heating,  whereas  the  second  term  is  regarded  with  the  work  pro¬ 
duced  by  the  electric  current  against  the  Seebeck  effect.  Invoking 
the  following  thermoelectric  relation  [15], 

zdT  =  T  doc  (6) 

and  substituting  Eqs.  (3)  and  (5)  into  Eq.  (1),  the  following  equation 
for  the  temperature  distribution  in  a  thermoelectric  material  can  be 
derived, 

V-(kVT)-rj-VT  +  p(J-j)=0  (7) 

where  the  first  term  refers  to  the  Fourier  (heat)  conduction,  the 
second  one  to  the  Thomson  (thermoelectric)  effect,  and  the  third 
one  to  the  Joule  heating.  Writing  Eq.  (4)  for  the  electric  current, 
recalling  that  V  j  =  0  at  steady-state  conditions,  one  can  derive 
the  following  expression  for  the  voltage  distribution  along  the 
domain, 

V  •  (yw)  +  V  •  (yaVT)  =  0  (8) 

where  y  =  1/p  is  the  electrical  conductivity.  The  first  term  stands  for 
the  electric  conduction,  whereas  the  second  one  refers  to  the 
distortion  on  the  electric  field  induced  by  the  thermoelectric 
effect.  Eqs.  (7)  and  (8)  rule  the  thermoelectric  phenomena,  being 
both  expressed  for  2-D  Cartesian  domain  as  follows: 


dx 


-  T 


f.m_  .  dT\ 
)xfa+Jydy) 


=  0 


0) 


—  (y— ) 

dx  y  dx ) 


—  (y— ) 

dy  dy) 


d  (  dT\  d  (  dT 

+  ¥xyadx)+dyyady 


=  0 


(10) 


where  jx  and  jy  are  the  x  and  y  components  of  the  current  density, 
respectively,  in  (A/m2).  Eqs.  (9)  and  (10)  require  two  boundary  con¬ 
ditions  each.  For  the  latter,  prescribed  inlet  (Vin)  and  outlet  ( Vout ) 
voltages  were  adopted.  In  addition,  bearing  in  mind  that  there  is 
no  electron  flux  through  subdomains  1  and  6,  dV/dy  =  0  boundary 
conditions  have  also  been  adopted.  In  case  of  Eq.  (9),  prescribed 
temperatures  were  used  for  both  hot  ( Th )  and  cold  ( Tc  =  Th-AT ) 
ends.  Zero  heat  flux  boundary  conditions  ( dT/dx  =  0)  were  addition¬ 
ally  employed  for  the  cell  symmetry.  Fig.  2  depicts  the  conditions 
used  for  each  boundary  of  the  physical  domain.  The  thermophysical 
properties  of  the  Bi2Te3-elements  were  calculated  from  2nd-order 
polynomial  fits  obtained  from  data  provided  by  Rowe  [16], 

p  =  2.3935  x  10“nT2  +  2.9771  x  10“8T  -  8.959  x  107  (11) 


k  =  3.682  x  lO'5!2  -2.372  x  lCT2r  +  5.388 


(12) 


a  =  -8.5952  x  10  1072  +  8.0546  x  KT7T  + 4.329  x  10“5  (13) 

where  T  is  in  (I<),  p  in  (Vm/A),  k  in  (W/mK),  and  a  in  (V/K).  The 
Thomson  coefficient,  t,  was  calculated  from  Eqs.  (6)  and  (13).  Both 
the  air  and  the  Al203-elements  were  assumed  to  be  perfect  electric 
insulators,  with  thermal  conductivities  (at  300  K)  of  30  and 
0.026  W/mK,  respectively.  For  the  copper,  a  thermal  conductivity 
of  400  W/mK  and  an  electrical  conductivity  of  5.96  x  107A/Vm 
have  been  adopted.  The  heat  transfer  inside  the  air  cavity  has  been 
modelled  assuming  a  unitary  Nusselt  number,  so  that  the  effects  of 
free  convection  and  radiation  have  been  neglected. 


Fig.  2.  Physical  domain  with  boundary  conditions. 
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Fig.  3.  Typical  control  volume  of  the  discretized  domain. 


volumes.  A  snapshot  of  a  computational  domain  with  3120  control 
volumes  is  depicted  in  Fig.  4. 

Eqs.  (9)  and  (10)  were  discretized  using  a  2nd-order  central- 
differencing  scheme.  The  non-linear  terms  have  been  incorporated 
into  the  source  term.  The  resulting  algebraic  equations  for 
temperature  and  voltage  are  as  follows: 

ApTp  =  ATe  Te  +  AtwTw  +  ATnTN  +  A(  Ts  +  BT 
'  Al  =  keAy/Sxe 
Aw  —  kwAy/Sxw 
ATn  =  knAx/Syn 
'  AST  =  ksAx/dys 
ATP  =  ATe  +  ATW  +  ATn  +  A l 

„  BT  =  pP(j2x  +j2y)AxAy  -  t P(jx(re  -  Vw)Ay  +jy(Tn  -  Vs) Ax) 

(14) 


3.  Numerical  scheme 

Because  of  the  non-linearities,  Eqs.  (9)  and  (10)  must  be  solved 
iteratively  in  order  to  come  out  with  the  temperature  and  voltage 
distributions  along  the  domain.  The  computational  code  was  writ¬ 
ten  in-house  based  on  the  so-called  finite-volume  method  [17]  to 
provide  access  to  the  numerical  schemes  in  the  present  and  also 
in  further  researches.  The  method  consists  of  dividing  the  physical 
domain  into  non-overlapping  control  volumes  in  which  the 
transported  quantities  are  conserved.  The  centroid  of  each  control 
volume,  as  illustrated  in  Fig.  3,  corresponds  to  an  integration  point 
of  the  discretized  domain.  The  properties  (T,  V)  are  evaluated  at  the 
centroids,  whereas  the  fluxes  (q,  j)  are  evaluated  at  the  control  sur¬ 
faces.  A  non-uniform  Cartesian  mesh  was  generated  by  means  of 
the  equation  introduced  by  Wood  [18].  Mesh  independent  solu¬ 
tions  have  been  found  for  computational  grids  with  3000+  control 


AvpVP  =  Ave  VE  +  AvwVw  +AvnVN+  Avs  Vs  +  Bv 

'  Ave  =  yeAy/Sxe 

Al  =  ywAy/Sxw 

K  =  y„Ax/<Syn 

<  Avs  =ysAx/dys 

Avp=Ave  +Avw+Avn+Avs 

Bv  =  CLeAverE  +  aw Alrw  +  a  nAvnrN  +  a  sA(rs 

-  (a eAve  +  ccwAl  +  aX  +  aAV )  T"P 

where  the  superscripted  asterisk  stands  for  the  property  evaluated 
from  the  previous  iteration.  The  components  of  the  electron  flux  at 
the  control  surfaces  n,  s,  e,  and  w  are  calculated  as  follows: 


Fig.  4.  Non-uniform  grid  with  3120  control  volumes  (see  Table  2  for  the  summary  of  the  physical  subdomains). 
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Fig.  5.  Continuity  of  heat  and  current  fluxes  at  the  interface. 
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Jn 


Fig.  8.  Comparison  between  calculated  and  experimental  COP. 

i^[(VN-VP)  +  an(TN-TP)}  (16) 


js  =  -^-[(VP-Vs)  +  as(TP-Ts)] 


(17) 


Fig.  6.  Comparison  between  calculated  and  experimental  electric  current:  (a) 
AT  =  0  K,  (b)  AT  =  30  K  and  (c)  A T=  60  K. 


je  =  -T^[(VE-VP)  +  ae(TE-TP)]  (18) 

We 

Jw  =  -  It  ttVp  -  -  Tw)l  (19) 

('Jw 

Therefore,  the  x  and  y  components  of  the  electron  flux  are  cal¬ 
culated  from: 

:  _ie^xe  +  jwSxw 

Jx~  dxe  +  dxw  1  j 


:  _  Jn^yn  Jsfos 

Jy~  Syn  +  3ys 


(21) 


The  temperature  and  voltage  at  the  interfaces  between  different 
materials  are  calculated  in  order  to  guarantee  the  continuity  of  the 
electric  current  and  heat  fluxes  (see  Fig.  5)  as  follows: 
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Fig.  9.  Temperature  distribution  along  the  domain  for  AV  =  16  V  and  different  AT  (1  and  5  -  A1203  wafer,  2  and  4  -  Cu  contacts  ,  and  3  -  thermoelectric  element). 
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Te  = 


(1  -fe)TP+feTE  + 


(alj)£  -  (ocTj)p 

2  kp  ,  2  kE 
Axp  '  Axe 


thermoelectric  effect  is  on,  the  locus  of  the  maximum  temperature 
migrates  from  the  centre  to  the  bottom  inasmuch  AT  increases. 


=  (1  -fe)VP  +feVE  +  a£(rp-g-g(Wp)  (23) 

where /e  =  lA xE/Sxe.  The  sets  of  linear  equations  have  been  solved 
iteratively  through  the  so-called  tri-diagonal  matrix  algorithm 
(TDMA).  In  addition  to  the  boundary  conditions,  the  input  variables 
are  the  thermophysical  properties  of  the  composing  materials,  and 
the  cell  dimensions,  all  described  in  Table  2.  More  detailed  informa¬ 
tion  on  the  numerical  scheme  can  be  found  in  [19]. 

4.  Results  and  discussion 

4.1.  Model  validation 

The  code  predictions  were  validated  against  experimental  data 
obtained  from  the  manufacturer  of  a  particular  thermoelectric 
module  [13].  All  simulations  were  carried  out  for  Th  =  323  K,  but 
varying  the  AT  between  the  hot  and  the  cold  ends  from  0  to 
60  K,  and  the  AV  applied  to  the  whole  thermoelectric  cell.  Fig.  6 
shows  a  comparison  between  the  calculated  and  measured  electric 
current,  where  one  can  see  the  maximum  difference  achieved  (for 
AT  =  0  K  and  AV  =  16  V)  was  below  the  10%  threshold,  while  the 
average  difference  lies  around  4%.  In  all  cases,  it  can  be  seen  the 
model  is  able  to  reproduce  the  experimental  trends  satisfactorily. 
One  should  also  note  that  Fig.  6  illustrates  the  behaviour  described 
by  Eq.  (4),  where  a  decreasing  AT  leads  to  a  current  increase  for  a 
given  A V  because  of  the  Seebeck  effect. 

Additional  validation  parameters  are  the  cooling  capacity  Qc 
and  the  COP,  calculated  respectively  as  follows: 

(24) 


4.3.  Thermoelectric  properties 


The  sensitivity  analysis  was  carried  out  considering  as 
response  variables  the  cooling  capacity  and  the  COP,  whereas  the 
thermophysical  properties  (i.e.  thermal  conductivity,  Seebeck  coef¬ 
ficient,  and  electric  conductivity)  have  been  taken  as  independent 
parameters.  A  2-level,  3-factor  factorial  design  was  then  structured 
totalizing  8  runs.  The  levels  span  ±5%  taking  the  figures  provided 
by  Eqs.  (11 )— (13)  as  reference.  The  simulation  runs  were  carried 
out  for  AT  =  0  K  and  AV=  16  V.  The  regression  model  adopted  in 
this  work  is  as  follows: 


W  —  Aq  H-  A\  h  2,2 dt  hy  T-  A4kdc  A$lcy  -|-  A^ocy  A^icoty  (27) 

where  7  =  1/p  is  the  electric  conductivity,  W  is  the  dimensionless 
response  variable,  A  are  the  coefficients  calculated  from  the  least- 
squares  method,  and  cp  are  the  dimensionless  values  of  cp ,  calcu¬ 
lated  from: 


q>  =  2-^ — 

*Pmax  —  *Pmin 


(28) 


Fig.  10  shows  the  cooling  capacity  is  mainly  affected  by  the 
electric  conductivity  and  the  Seebeck  coefficient,  and  marginally 
affected  by  the  thermal  conductivity,  which  play  a  negative  role 
on  the  cooling  capacity.  Higher  order  interactions  have  not  played 
any  material  roles  on  the  cooling  capacity.  Fig.  10  also  shows  the 
effects  of  the  thermophysical  properties  on  the  COP,  where  one 
can  see  the  Seebeck  coefficient  plays  a  dominant  role,  followed 
by  the  thermal  and  electric  conductivities,  which  played  marginal 
effects.  These  results  are  confirmed  by  the  definition  of  the  figure- 
of-merit  of  the  thermoelectric  material,  Z=  a2y/k,  which  is  straight¬ 
forwardly  related  to  the  COP  15]. 


4.4.  Aspect  ratio 


COP  =  .  ^  .  (25) 

Qi,-Qc 

where  Qh  is  the  heat  released  at  the  hot  end,  calculated  as  follows: 

Qh=NLz^(fc7A^rlAx)  (26) 

where  Lz  is  the  cell  dimension  in  the  z-direction  (perpendicular  to 
the  paper  sheet),  N  is  the  number  of  thermoelectric  couples  in  the 
cell,  and  n  and  m  are  the  number  of  control  volumes  in  the  x  and  y 
directions,  respectively.  Fig.  7  shows  the  calculated  and  the  experi¬ 
mental  cooling  capacities  agreed  to  with  errors  within  the  10% 
threshold.  The  higher  difference  is  observed  for  low  voltages  and 
AT  =  0  K.  A  similar  behaviour  is  observed  in  Fig.  8  for  the  COP.  In 
all  cases,  the  experimental  trends  are  well  reproduced  by  the  model. 

4.2.  Temperature  distribution 

Fig.  9  explores  the  temperature  distributions  obtained  for  four 
different  cases:  (a)  no  thermoelectric  effect  (Joule  heating  only) 
and  AT  =  0  K;  (b)  thermoelectric  effect  and  AT  =  0  K;  (c)  thermo¬ 
electric  effect  and  AT  =  30  K;  and  (d)  thermoelectric  effect  and 
AT  =  60  K.  In  all  cases,  AV=  16  V.  The  temperature  profiles  along 
the  A-A  cut  (at  x  =  1.1  mm)  are  also  depicted  in  Fig.  9.  For  case  (a), 
where  no  thermoelectric  effect  takes  place,  one  can  see  that  the 
Joule  heating  is  symmetrically  dissipated  by  Fourier  conduction  in 
such  a  way  the  maximum  temperature  takes  place  at  the  centre 
of  the  thermoelectric  elements.  This  is  so  as  AT  =  0  K.  In  cases  the 


To  assess  the  influence  of  the  geometry,  the  aspect  ratio  was 
varied  by  increasing  the  height  of  the  thermoelectric  cell,  Ly,  in 
two  fashions:  (a)  constrained  base  area  (i.e.  fixed  Lx,  see  Fig.  11a), 
and  (b)  constrained  volume  of  thermoelectric  material  (see 
Fig.  l  ib).  In  all  cases,  AT  =  0  K  and  the  voltage  was  varied  from 
14  to  20  V.  Fig.  12  shows  the  COP  is  weakly  affected  by  Ly.  Actually, 
a  slight  increase  can  be  observed.  This  is  so  as  the  cooling  capacity 
depletes  inasmuch  the  electric  current  decreases,  which  dimin¬ 
ishes  the  power  consumption  at  the  same  rate.  As  the  COP  stands 
for  the  ratio  between  the  cooling  capacity  and  the  power  consump¬ 
tion,  one  may  expect  the  COP-value  does  not  change  significantly 
from  one  case  to  the  other. 


Sensitivity 

Fig.  10.  Results  of  the  sensitivity  analysis:  cooling  capacity  and  COP. 
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Fig.  11.  Samples  of  geometries  analysed  in  case  of  (a)  constrained  base  area,  and  (b)  constrained  volume. 
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Fig.  12.  Influence  of  the  aspect  ratio  on  the  COP  in  case  of  constrained  area  (solid 
bullets)  and  constrained  volume  (open  bullets). 


5.  Summary  and  conclusions 

A  two-dimensional  steady-state  model  was  put  forward  to  eval¬ 
uate  the  thermodynamic  performance  of  thermoelectric  cells  in  the 
realm  of  the  irreversible  thermodynamics.  The  model  takes  into 
account  the  Fourier  conduction,  the  Joule  heating,  and  the 
Thomson  effect,  being  able  to  predict  the  cooling  capacity,  the 
power  consumption,  and  the  COP  in  case  of  prescribed  voltage 
supply,  and  prescribed  temperatures  at  the  hot  and  cold  ends. 
The  governing  equations  were  discretized  by  means  of  the 
finite-volume  method  using  a  central-differencing  scheme.  The 
non-linearities  typical  of  the  thermoelectric  phenomena  were 
embedded  into  the  source  term,  and  the  resulting  sets  of  algebraic 
equations  were  solved  iteratively  by  the  TDMA  algorithm. 

The  tailor-made  model  was  coded  in-house  and  its  predictions 
for  electric  current,  cooling  capacity  and  COP  were  compared 
against  experimental  data  obtained  from  the  manufacturer  of  a 
particular  thermoelectric  cell.  It  was  observed  that  the  numerical 
predictions  and  experimental  data  not  only  agreed  to  within  8.8% 
thresholds  (electric  current),  10%  thresholds  (cooling  capacity), 
and  7.3%  thresholds  (COP),  but  also  that  the  model  is  able  to  follow 
the  experimental  trends  very  closely.  The  model  can  be  used  either 
to  access  the  performance  of  an  existing  cell  for  given  working 


conditions  (AT,  AV  or  current)  or  to  design  a  new  one  (geometry, 
selection  of  materials)  subjected  to  working  constraints  (cooling 
capacity,  COP). 

The  influence  of  the  thermophysical  properties  on  the  response 
variables  (cooling  capacity  and  COP)  was  assessed  by  means  of  a  23 
factorial  design,  which  has  pointed  out  that  the  Seebeck  coefficient 
and  the  thermal  conductivity  play  major  roles  on  the  cooling 
capacity,  whereas  the  COP  is  more  sensible  to  the  Seebeck  coeffi¬ 
cient.  The  influence  of  the  geometry  was  also  assessed  by  varying 
the  aspect  ratio  according  to  two  different  ways:  constrained  base 
area  and  constrained  volume.  It  was  found  that  both  the  cooling 
capacity  and  the  power  consumption  vary  at  the  same  rate,  in  such 
a  way  the  COP  has  showed  a  comparable  behaviour  for  constrained 
base  area  and  constrained  volume  of  thermoelectric  material. 
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